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This lecture contains a brief introduction to HERA and deep inelastic scattering (DIS), 
before going on to highlight some of the measurements of the hadronic final state in DIS 
— «, ' performed by the HI and ZEUS collaborations. 

ix. 



1 The HERA Accelerator and Detectors 



The HERA accelerator, located at DESY in Hamburg, is an electron-proton collider. It is 
6.3 km in circumference and collides positrons (or electrons) at pa 27 GeV with protons at 
820 GeV. There are four interaction regions: two containing general purpose, hermetic detectors 
(HI and ZEUS); another experiment (HERMES) investigating the spin distributions of the 
quarks in protons and neutrons; and another (HERA-B) planning to measure CP violation in 
the B-system. The HI and ZEUS detectors took first data in 1992. 

The ZEUS detector is shown in figure |l| The asymmetric design of the detector reflects the 
proton energy being significantly higher than that of the electron beam. 



Overview of the ZEUS Detector 
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Figure 1: Cross sectional view of the ZEUS detector 

The tracking system consists of a vertex detector (VXD) [|TJ and a central tracking cham- 
ber (CTD) || enclosed in a 1.43 T solenoidal magnetic field. Immediately surrounding the 
beampipe is the VXD which consists of 120 radial cells, each with 12 sense wires. The CTD, 
which encloses the VXD, is a drift chamber consisting of 72 cylindrical layers, arranged in 9 
superlayers. Superlayers with wires parallel to the beam axis alternate with those inclined at 
a small angle to give a stereo view. A forward tracking detector is employed in the forward 
region to detect tracks in the proton direction and consists of three 12-layer planar drift cham- 
bers sandwiched with pairs of transition radiation detectors. In the rear direction there is an 
additional 12-layer planar drift chamber known as the rear tracking detector (RTD). 

Outside the solenoid is the uranium-scintillator calorimeter (CAL) which is divided into 
three parts: forward, barrel and rear covering the polar regions 2.6° to 36.7°, 36.7° to 129.1° and 
129.1° to 176.2°, respectively. The CAL covers 99.7% of the solid angle, with holes of 20 x 20 
cm 2 in the centres of the forward and rear calorimeters to accommodate the HERA beam pipe. 
Each of the calorimeter parts is subdivided into towers which are segmented longitudinally into 
electromagnetic (EMC) and hadronic (HAC) sections. The small angle rear tracking detector 
(SRTD) Q, which is attached to the front face of the rear calorimeter, measures the impact 
point of charged particles at small angles with respect to the positron beam direction. 

The iron return yoke for the magnet is instrumented with proportional counters. This 
backing calorimeter (BAC) measures any hadronic energy which 'leaks out' out of the main 
calorimeter. Beyond that and in the forward direction there are further detectors for muon 
detection. 

Downstream of the main detector in the proton direction, six measuring stations are installed 
in the proton ring for detecting forward scattered protons. Beyond the final station, further 
downstream, is a forward neutron calorimeter. In the electron direction, two lead scintillator 



calorimeters placed —35 m and —107 m from the interaction point measure the luminosity and 
tag events with a small momentum transfer []5| . 

A fuller description of the ZEUS detector can be found in reference @. The HI detector is 
of a very similar layout as ZEUS and a description can be found in reference 0. 

2 DIS Kinematics 

The event kinematics of deep inelastic scattering, DIS, are determined by the negative square 
of the four-momentum transfer at the positron vertex, Q 2 = —q 2 , and the Bjorken scaling 
variable, x = Q 2 /2P ■ q, where P is the four-momentum of the proton. In the Quark Parton 
Model (QPM), the interacting quark from the proton carries the four-momentum xP. The 
variable y, the fractional energy transfer to the proton in its rest frame, is related to x and Q 2 
by y ~ Q 2 /xs, where ^/s is the positron-proton centre of mass energy. Because the HI and 
ZEUS detectors are almost hermetic the kinematic variables x and Q 2 can be reconstructed in 
a variety of ways using combinations of electron and hadronic system energies and angles || . 
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Figure 2: (a) QPM (b) QCDC and (c) BGF diagrams 

In QPM there is a 1+1 parton configuration, fig. ^|a, which consists of a single struck quark 
and the proton remnant, denoted by "+1". At HERA energies there are significant higher- 
order Quantum Chromodynamic (QCD) corrections: to leading order in the strong coupling 
constant, a s , these are QCD-Compton scattering (QCDC), where a gluon is radiated by the 
scattered quark and Boson-Gluon- Fusion (BGF), where the virtual boson and a gluon fuse to 
form a quark-antiquark pair. Both processes have 2+1 partons in the final state, as shown in 
fig. 0. There also exists calculations for the higher, next-to-leading (NLO) processes. 

Perturbative QCD does not predict the absolute value of the parton densities within the 
proton but determines how they vary from a given input. For a given initial distribution at a 
particular scale Altarelli-Parisi (DGLAP) evolution || enables the distributions at higher Q 2 
to be determined. DGLAP evolution resums the leading log(Q 2 ) contributions associated with 
a chain of gluon emissions. At large enough electron-proton centre-of-mass energies there is a 
second large variable 1/x and, therefore, it is also necessary to resum the log(l/x) contributions. 



This is acheived by using the BFKL equation [10 



3 Jet Physics 

To relate the hadronic final state to the underlying hard partonic behaviour it is generally 
necessary to apply a jet algorithm. The JADE algorithm jTI]] has been used in the following 
analyses as it was, at the time, the only algorithm which allowed comparison to the NLO 
calculations (PRO JET @ and DIS JET 0). The JADE algorithm is a cluster algorithm 



based on the scaled invariant mass-squared 



JADE _ 2EjEj(l - cos%) 

for any two objects i and j assuming that these objects are massless. W 2 is the squared invariant 
mass of the hadronic final state and 9ij is the angle between the two objects of energies E{ and 
Ej. The minimum yij of all possible combinations is found. If the value of this minimum 
Uij is less than the variable cut-off parameter y cut , the two objects i and j are merged into 
a new object by adding their four-momenta and the process is repeated until all y^ > y cut . 
The surviving objects are called jets which represent the underlying partonic structure that is 
dependent on a s . 




Figure 3: Jet production rates Rj as a function of the jet resolution parameter y cut for Q in the 
range (a) 120 < Q 2 < 240 GeV 2 , (b) 240 < Q 2 < 720 GeV 2 , (c) 720 < Q 2 < 3600 GeV 2 , and 
(d) 120 < Q 2 < 3600 GeV 2 . Only statistical errors are shown. Two NLO QCD calculations, 
DISJET and PRO JET, each with the value of Aj^ obtained from the fit at y cut =0.02, are also 
shown. 

Figures §a-d show the ZEUS jet rates using data taken in 1994, R\+\, R2+1 and as a 
function of y cut for data compared with the DISJET and PROJET NLO QCD calculations for 
three Q 2 intervals 120 < Q 2 < 240 GeV 2 , 240 < Q 2 < 720 GeV 2 , 720 < Q 2 < 3600 GeV 2 , and 
the combined region 120 < Q 2 < 3600 GeV 2 . There is good agreement between the corrected 
1 + 1 and 2 + 1 jet rates and the NLO QCD calculation over most of the range in y cut shown. 
Both programs agree well in their prediction of the jet-rate dependence as a function of y cut . 
The values of ac s (Q) extracted by the HI ]TJ| and ZEUS |T5j collaboration as a function of 



Q are shown in Fig. f|. The value of a s was determined by varying the A scale parameter in 
the QCD calculation until the best fit to the ratio R2+1 was obtained at a particular value of 
y cut . The measured a s decreases with increasing Q, consistent with the running of the strong 
coupling constant, with Q 2 taken as the scale. In addition the figure shows the curves for 
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Figure 4: Left: Values and total error ofa s (Mz) from various processes. The solid line indicates 
the world average and the band its total error. Right: a s (Q) from HERA (open symbols) 
and other processes with increasing Q (closed circles): T T and Ohad/otot? event shapes and 
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100, 200, and 300 MeV. An extrapolation to a s {M z ) yields: 



HI 93 : a s (M z ) = 0.123 ± 0.012(stat) ± 0.013(syst.) 
ZEUS 94 : a s (M z ) = 0.117 ± 0.005(stat)+J;J^(exp.) ± 0.007(th.) 

which are consistent with other values obtained from a large variety of different processes as 
shown Fig. |] (see [16] for references). Even with the current statistics the HERA measurements 
are already competitive with those made elsewhere. 

Recently two new, more flexible NLO calculations (MEPJET fL7| and DISENT |J|]) have 
become available allowing the experiments to analyze the data using any particular jet algo- 
rithm. The algorithm |l9j is particularly suited for DIS as it allows factorization between 

The ZEUS collaboration has reanalyzed |2T 



the beam fragmentation and the hard process 
their 1994 data using this algorithm. The preliminary values of a s (Q) obtained in the three 
bins of Q are shown (with statistical errors only) in Fig. |3] and are consistent with the results 
obtained with the JADE algorithm. 



4 Event Shapes 

A natural frame in which to study the dynamics of the hadronic final state in DIS is the Breit 
frame p2]| . In this frame the exchanged virtual boson is purely space-like with 3- momentum 
q = (0,0,— Q), the incident quark carries momentum Q/2 in the positive Z direction, and 
the outgoing struck quark carries Q/2 in the negative Z direction. A final state particle has a 
4-momentum p B in this frame, and is assigned to the current region if is negative, and to the 
target frame if pf is positive. The advantage of this frame lies in the maximal separation of the 
outgoing parton from radiation associated with the incoming parton and the proton remnant, 
thus providing the optimal environment for the study of the fragmentation of the outgoing 
parton. 

Event shape variables have been investigated in e + e~ experiments and used to extract the 



strong coupling constant a s (Mz) independent of any jet algorithm, see eg ref . [23]. HI have re- 
cently performed a similar analysis in deep inelastic scattering in the current fragmentation 
region of the Breit frame. 



The event shape dependence on Q (or energy dependence) can be due to the logarithmic 
change of the strong coupling constant a s (Q) oc 1/lnQ, and/or power corrections (hadronisa- 
tion effects) which are expected to behave like 1/Q. Recent theoretical developments suggest 
that 1/Q corrections are not necessarily related to hadronisation, but may instead be a uni- 
versal soft gluon phenomenon associated with the behaviour of the running coupling at small 
momentum scales [^5|. |2f| . 

HI have analysed a number of infrared safe (ie independent of the number of partons 
produced) event shape variables. Their definitions are given below, where the sums extend over 
all hadrons h (being a calorimetric cluster in the detector or a parton in the QCD calculations) 
with four-momentum pf = {E^, pf } The current hemisphere axis n = {0, 0, —1} coincides 
with the virtual boson direction. 
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T z = h 5 — = — — : — ^— n = hemisphere axis 



• Jet Broadening B c 
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A common characteristic of the mean event shape values (1 — T c ), (1 — T z ), ( B c ) and (p c ) 
is the fact that they exhibit a clear decrease with rising Q, fig. [5]. This is due to fact that the 
energy flow becomes more collimated along the event shape axis as Q increases, a phenomenon 
also observed in e + e~ annihilation experiments. 

HI showed by fitting to the data in fig. |5| all the event shape variables can be well described by 
just the first order power corrections oc 1/Q, without the need for any higher order corrections. 
The second order perturbative QCD parton predictions are also shown and their discrepancies 
with the data show that the power corrections are substantial at low values of Q, but become 
less important with increasing energy. 

The analysis of the event shapes give results consistent with each other for ato, the power 
correction parameter thus supporting the prediction of universality p5| , and also gives consis- 
tent values of a s (M z ). The results of the fit are a = 0.491 ± 0.003 (exp) 1q!o42 (theory) for the 
power correction parameter and a s (Mz) — 0.118 ± 0.001 (exp) 1q!oo6 (theory) for the strong 
coupling constant in the MS scheme. These values are compatible with those extracted by e + e~ 
experiments |?7| 
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Figure 5: Mean event shape variables as a function of Q for a) (1 — T C ), b) (l — T z )/2, c) (B c ), 
and d) (p c ) . HI DIS ep data (•, errors include statistics and systematics) are compared with 
Q CD fits ( — ) and second order QCD calculations (■ ■ ■) 



5 Fragmentation Functions 

Fragmentation functions represent the probability for a parton to fragment into a particular 
hadron carrying a certain fraction of the parton's energy. Fragmentation functions incorporate 
the long distance, non-perturbative physics of the hadronization process in which the observed 
hadrons are formed from final state partons of the hard scattering process and, like structure 
functions, cannot be calculated in perturbative QCD, but can be evolved from a starting dis- 
tribution at a defined energy scale. If the fragmentation functions are combined with the cross 
sections for the inclusive production of each parton type in the given physical process, pre- 
dictions can be made for the scaled momentum, x p , spectra of final state hadrons. Small x p 
fragmentation is significantly affected by the coherence (destructive interference) of soft glu- 
ons [2£], whilst scaling violation of the fragmentation function at large x p allows a measurement 

ofa s [0- 

In e + e~ annihilation the two quarks are produced with equal and opposite momenta, ±a/s/2. 
This can be compared with a quark struck from within the proton with outgoing momentum 
—Q/2 in the Breit frame. In the direction of the struck quark (the current fragmentation 
region) the particle momentum spectra, x p = 2p B /Q, are expected to have a dependence on Q 
similar to those observed in e + e~ annihilation [H], [!2j at energy yfs = Q. 

The inclusive charged particle distributions |33], |34j| , (l/a to t)da/dx p , are shown in figure || 
plotted in bins of fixed x p as a function of Q 2 . For Q 2 > 80 GeV 2 the distributions rise with 
Q 2 at low x p and fall-off at high x p and high Q 2 . By measuring the amount of scaling violation 
one can ultimately measure the amount of parton radiation and thus determine a s . Below 
Q 2 = 80 GeV 2 the fall off is due to depopulation of the current region. 

The results can be compared to the next-to-leading order (NLO) QCD calculations, as 
implemented in CYCLOPS |j3E| , of the charged particle inclusive distributions in the restricted 
region Q 2 > 80 GeV 2 and x p > 0.1, where the theoretical uncertainties are small, unaffected by 
the hadron mass effects which are not included in the fragmentation function. This comparison 
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Figure 6: The inclusive charged particle distribution, l/o~ to t da/dx p , in the current fragmenta- 



tion region of the Breit frame compared to the NLO calculation, CYCLOPS MW. 



is shown in figure |]. The NLO calculation combines a full next-to-leading order matrix element 
with the MRSA' parton densities (with a Aqcd = 230 MeV) and NLO fragmentation functions 
derived by Binnewies et al. from fits to e + e~ data ||36|| . The data and the NLO calculations 
are in good agreement, supporting the idea of universality of quark fragmentation. 

The peak position of the £ = \n(l/x p ) distributions, £ pea k, was evaluated. Figure |7| shows 
the distribution of £ pca k a s a function of Q for the HERA data [33, 37, [38|] and of y/s for the 
e + e~ data. Over the range shown the peak moves from ~ 1.5 to 3.0, equivalent to the position 
of the maximum of the corresponding momentum spectrum increasing from ~ 400 to 900 MeV. 



The HERA data points are consistent with those from TASSO |39[ data and a clear agreement 
in the rate of growth of the HERA points with the e + e" 



data |39|, [4(| is observed. 
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Figure 7: £ pea k a function ofQ. The HERA data are compared to results from OPAL, TASSO 
and TOPAZ. A straight line fit of the form £ pea k = b ln(Q) + c to the ZEUS £ pea k values is 
indicated as well as the line corresponding to b = 1, discussed in the text. 



The increase of £ pea k can be approximated phenomenologically by the straight line fit £ pea k = 
b ln(Q) + c also shown in figure 0. Also shown is the statistical fit to the data when 6 = 1 which 



would be the case if the QCD cascade was of an incoherent nature, dominated by cylindrical 
phase space. (A discussion of phase space effects is given in ]4TJ| .) In such a case, the logarithmic 
particle momentum spectrum would be peaked at a constant value of momentum, independent 
of Q. The observed gradient is clearly inconsistent with 6=1 and therefore inconsistent with 
cylindrical phase space thus supporting the coherent nature of gluon radiation. 



6 BFKL versus DGLAP Evolution 

In the DGLAP parton evolution scheme @ the parton cascade follows a strong ordering in 
transverse momentum k\ n 3> k\ n _ x ^> ... ^> while there is only a soft (kinematical) 
ordering for the fractional momentum x n < x n -\ < ... < x\. However for low-x at HERA the 
BFKL scheme [1111 could well be the dominant scheme. In this scheme the cascade follows a 



strong ordering in fractional momentum x n <C x n -\ <C ... <C x%, while there is no ordering in 
transverse momentum. 

BFKL evolution can be enhanced by studying DIS events which contain an identified jet of 
longitudinal momentum fraction Xj et = p z (j)/E prot on (in the proton direction) which is large 
compared to Bjorken x [42[. By tagging a forward jet with pr{j) — Q this allows minimal phase 
space for DGLAP evolution while the condition Xj e t 3> x leaves BFKL evolution active. This 
leads to the forward jet production cross section in BFKL dynamics being larger than that of 
the 0(a%) QCD calculation with DGLAP evolution J43[. 

In Fig. H recent data from HI f|4| and ZEUS |45| are compared with BFKL predictions |46 



and fixed order QCD predictions as calculated with the MEPJET [TjJ program at NLO. The 
conditions pr{j) — Q and Xj e t 3> x are satisfied in the two experiments by slightly different 
selection cuts. HI selects events with a forward jet of pr(j) > 3.5 GeV (in the angular region 
7° < 0{j) < 20°) with 

0.5 < p T (j) 2 /Q 2 < 2 , x jet ~ E jet /E protm > 0.035 ; (1) 

while ZEUS triggers on somewhat harder jets of pr{j) > 5 GeV and < 2.4 with 

•Ejet Pz 

{j)/E proton > 0.035 . (2) 

Fig. shows that both experiments observe substantially more forward jet events than 
expected from NLO QCD. A very rough estimate of the uncertainty of the NLO calculation 
is provided by the two dotted lines, which correspond to variations by a factor 10 of the 
renormalisation and factorisation scales [<? R and [j? F . A recent BFKL calculation (dashed lines) 
gives a better agreement with the data. The overall normalisation in this calculation is uncertain 
and the agreement may be fortuitous, indeed it should also be noted that both experiments 
observe more centrally produced dijet events than predicted by the NLO QCD calculations. 
Further careful investigation is necessary before claiming that BFKL is the mechanism for this 
enhanced forward jet production. 



7 Conclusions 

To understand the underlying QCD processes in DIS it is necessary to study the hadronic final 
state. At the current level of understanding, QCD works well and describes the HERA data. As 
the precision of the HERA data improves and the NLO QCD calculations become available the 
framework of QCD is being tested more thoroughly. As yet it is not possible to say conclusively 
whether the effects of BFKL dynamics are being observed in the HERA data, much theoretical 




Figure 8: Forward jet cross section at HERA as a function of Bjorken x within (a) the HI 
and (b) the ZEUS j^/ acceptance cuts. The BFKL result of Bartels et al. fiffij is shown as the 
dashed line. The solid and dotted line give the NLO MEPJET result and a measure for the 
uncertainity of NLO prediction through changes in the choice of scale. 



and experimental work is in progress to define and measure variables that will allow a definitive 
statement. 
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